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Abstract 

Renormalization factors for focal vector and axial vector currents for the 
Wilson quark action are perturbatively calculated to one loop order includ- 
ing finite quark masses from the ratio of the on-shell quark matrix elements 
in the Feynman gauge defined on the lattice and in the continuum. For large 
quark masses of order unity in lattice units, we find that finite quark mass 
effects are quite large: one-loop coefficients of the renormalization factors 
differ by 100% compared to those in the massless limit. 
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1 Introduction 



Calculation of weak matrix elements for heavy-light and heavy-heavy mesons 
represents a subject of great interest in lattice QCD, which is in principle capable 
of a precise determination of the matrix elements from the first principles. The 
main source of systematic errors standing in the way to this goal is large m^a 
corrections for heavy quark mass rriq in units of the lattice spacing a. In cur- 
rent numerical simulations using the Wilson quark action the magnitude of the 
6-quark mass in lattice units is of order m^a ^1 — 2. To control the large rriqa 
error lattice studies have to address two problems. One is improvement of the 
quark action to reduce cut-off effects following either the Symanzik approach 
or Wilson's renormalization group approach]^. Another problem is a precise cal- 
culation of renormalization factors which relate operators on the lattice to those 
in the continuum for massive quarks. This calculation may be pursued either by 
perturbative methods]^ or by non-perturbative one|^. An improvement of the 
perturbative calculation including the finite rriqa corrections is the subject of this 
article. 

Renormalization factors connecting the lattice operator to the continuum one 
consists of the wave-function part and the vertex part. For the wave-function part 
Kronfeld and Mackenzie argued that the tree level normalization of the on-shell 
wave-function for the Wilson quark action suffers from large rriqa corrections for 
heavy quark[^. At the one-loop level, the question of rriqa corrections to the quark 
self-energy has been addressed by Kronfeld and Mertens for the Wilson quark 
action in Ref . . For the vertex part there are no rUga corrections at the tree level 
on the lattice. At the one-loop level the vertex part has been calculated only in 
the massless limit so far. For this reason analyses of weak matrix elements so far 
have to employ the one-loop expression in the massless limit for the renormalization 
factors combined with the Kronfeld-Mackenzie normalization at the tree level even 
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for heavy quark masses. For light quark masses m^a <^ 1 we have no doubt that 
the massless expressions are a good approximation, while for the case of heavy 
quark masses m^a 1 we may naturally suspect that corrections depending on 
TUqa is large even at the one-loop level. 

The purpose of this paper is to complete the one-loop calculation of the renor- 
malization factors for the vector and axial vector currents for finite quark masses 
with the Wilson quark action, and to use the results to investigate the magnitude 
of TUqa corrections at the one-loop level. We calculate the renormalization factors 
of the bilinear operators from the ratio of the on-shell quark matrix elements in the 
Feynman gauge defined in the lattice regularization scheme and in the continuum 
for various combinations of the two external quark masses. For the continuum 
regularization scheme we employ the naive dimensional regularization(NDR) with 
the MS subtraction. In order to regularize infrared(IR) divergences which are 
generated in one-loop contributions for the on-shell wave-function renormalization 
factor and the vertex corrections, we supply a fictitious mass A to the gluon prop- 
agator both for the lattice scheme and the continuum one. For the lattice scheme 
we extract the one-loop terms independent of A following the method used for the 
massless case in Ref. J^. 

This paper is organized as follows. In Sec. 2 we give the lattice and continuum 
Feynman rules, and describe the strategy for one-loop calculation of renormal- 
ization factors of the bilinear operators for finite quark masses. In Sec. 3 we 
demonstrate the technique to extract the one-loop terms independent of A for 
the lattice on-shell wave-function renormalization factor, and present one-loop re- 
sults for the relation between the lattice and continuum on-shell wave-function 
renormalization factors. Results for one-loop relations between the lattice vertex 
functions and the continuum ones are given in Sec. 4. In Sec. 5 we present results 
for the renormalization factors of the vector and axial vector currents and discuss 
the magnitude of rUqa corrections at the one-loop level. In Sec. 6 we compare our 
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results for the renormalization factor of the heavy-hght axial vector current with 
the previous static results in the heavy quark mass region toward the static limit. 
Our conclusions are summarized in Sec. 7. 

Throughout this paper we use the same notation for quantities defined on the 
lattice and their counterparts in the continuum. However, in case of any possibility 
of confusion, we shall make a clear distinction between them. 

2 Formalism 
2.1 Feynman rules 

The partition function of the lattice theory defined on a four-dimensional Eu- 
clidean space-time lattice with lattice spacing a is given by 

I U,,^VU^{x)U,jVjj;f{x)Vjj;f{x)e^p [So - Sw] , (1) 

where sites are labeled by x = (7110,^20,^30,^40) with ni,---,n4 integers. We 
take the standard Wilson gauge action for SU(3) gauge hnk variables Un{x) given 

by 

Sg = -\T. R-eTr [u^{x)U,{x + fi)Ul{x + i))Ul{x)] (2) 

with the bare coupling constant g. For the quark fields and V'/(a;) the Wilson 

quark action is given by 

Sw^a^Y.\ (^/ + t) ^/(^)^/(^) - ^ E [V'/l^^ + /i) (r + ^^)Ul{x)iljf{x) 

xj O / ZO ^ 

+V^/(a;)(r-7^)[/^(x)^/(a; + A)]}, (3) 

where /i is a vector with length a pointing along the /^-direction, ruf is the bare 
quark mass for each flavor / and r is the Wilson parameter. Color and spin 
indices are suppressed. We define the Euclidean gamma matrices in terms of the 
usual Minkowski matrices in the Bjorken-Drell convention according to 70 = "Jbd^ 
Ij = -^BD^ 75 = 7bd; they obey {7^, 7^} = 25^^ and 7J = 7^. 
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Gauge link variables are elements of SU(3) group in the fundamental represen- 
tation. They can be written in the form 



where Ta {A = 1, ■ ■ ■ , 8) are the generators of SU{3) group in the fundamental 
representation, which are normalized by Tt{TaTb) = 5ab/'^^ and G'^(x) are the 
gluon fields. In order to derive the Feynman rules we expand the link variable 
Ufj^{x) in terms of the coupling constant g. Higher order terms of iagJ^ATAG^^x) 
in the expansion of U^{x) yield tadpole graphs, if powers of G^{x) are contracted 
with each other. These tadpole contributions are suppressed only by powers of 
because of cancellations between powers of a from the expansion and ultraviolet 
divergences. As a result, coefficients in the perturbative expansion in g"^ are large, 
and lattice perturbation theory does not converge well. To avoid this problem we 
isolate the tadpole contributions as an overall constant Uq for the expansion of the 
link variable f/^(x) in terms of g: 



where uq is the expectation value of the link operator in the gauge employed for 
perturbative calculations 0. 

From now on throughout this paper we use lattice units for expressing physical 
quantities and suppress the lattice spacing a unless necessary. 

We summarize the lattice Feynman rules as follows. The gluon propagator in 
the Feynman gauge with momentum k is given by 



f/^(x)=exp iagY^TAG^^ix) 



(4) 



A 



U^{x)=uo l + iagY,TAG^{x) + ©(a^) 



(5) 



A 



^AB^^yUlDik, A), 



(6) 



where 



-1 



4^sin2(A;,/2) + A2. 



(7) 



a 
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Here, we give the gluon a small mass A, where eventually A ^ 0, to regularize pos- 
sible IR divergences in one-loop diagrams. The quark propagator with momentum 
k takes the form 



S ^{k,mu,r) = uo 



(8) 



with niu — {m + 4r — 4ruo)/uo. The one-gluon vertex with incoming quark mo- 
mentum p and outgoing momentum q has the following expression, 

-gTAUov^{p/2 + q/2,r), (9) 

where 

Vi^{k, r) = i7^cos(/c^) + rsin(/c^) (10) 

with no sum over At the one-loop level the two-gluon vertex appears only 
through gluon tadpole diagrams whose contributions are included in uq. 

The corresponding continuum Feynman rules are as follows: the gluon propa- 
gator in the Feynman gauge is 5AB5nvD{k, A) with D{k, \)~^ — k"^ -\- A^, the quark 
propagator is given by S~^{k,m) = {il/t + m), where I/: denotes J2alaka, and the 
quark-gluon vertex is —gT^v^ with — i^^. 

2.2 Procedure of calculation 

In the lattice regularization scheme ultraviolet divergences of composite opera- 
tors are regularized by the cutoff a~^, while in the NDR scheme in the continuum 
this is achieved by a reduction of the space-time dimension from four. Opera- 
tors defined in each regularization scheme can be related by renormalization fac- 
tors which are expected to converge for the perturbative expansion in terms of 
as — I (47r) due to the asymptotic freedom of the theory. For the vector and 
axial vector currents the relation takes the form, 

(^2(a:)7M^i(a:))^°''* = Zv^(;^2{x)-ii,'^\{x)^''^^ , (11) 
(^2 (a;) 7^75-01 (a;)) ^"""^ = (2;) 7^75^1(2;) j^''" , (12) 
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with Zy^ and Za^ (/x = 1, ■ ■ ■ , 4) renormalization factors. 

A possible way to perturbatively determine Zy^ and Za^ is to calculate the 

ratio of the matrix elements for the lattice and continuum regularization schemes 

employing external on-shell quark or anti-quark states in the Feynman gauge. 

For the vector and axial vector matrix elements we further need to specify space 

momenta of external quark or anti-quark states since the matrix elements generally 

depend on them. We take the natural choice of zero spatial momentum in this 

article. The renormalization constants are then calculated from 

7 _ (g2|V^2(x)7.^i(a;)|gi)-^ ^- - i 9 (^'^^ 

- (g2|^2(x)7.V^i(x)|gi)'- ^^"1'^'^^' ^^^^ 

^ (g2|^2(x)74V^i(x)|gi)^°'^^ ^^^^ 
(g2|V'2(a;)74^i(a;)|gi)^''" ' 

7 _ (g2|^2(x)7»75V'i(x)|gi)^°°' , . , . 

W2 1^2 (a;)7j75 V"! (a;) I gi) 

^ (g2|^2(3:)7475^l(a^)|gi)'°°' ^^g^ 

(g2|V^2(a;)7475V'i(a;)|gi)^''" ' 
where q and q stand for quark and anti-quark respectively. The choice of a quark 

state (^2! or an anti-quark state (52! for the external state is made to ensure a non- 
zero value of the matrix element for each operator. We should note that because 
of violation of space-time permutation symmetry in our choice of momenta for 
external quark or anti-quark states Zy. 7^ Zy^ and Za^ 7^ Za^ are expected due to 
possible iriqa corrections except in the limit of a ^ 0. 

At the tree level on-shell wave-function renormalization factor for the Wilson 
quark action is shifted from unity by finite m^a corrections]^ 
7(0) _ 7(0) _ 7(0) _ 7(0) 



= UoVcosh(E{°^) + rsinh(E{°Vcosh(4°^) + rsinh(4°^), (17) 

where E^^^ denotes the pole mass at the tree level in common between the Wilson 
quark action and the continuum one. The superscript (i) refers to the i-th loop 
level. Up to the one-loop level renormalization factors are written as 

Zy = 4°)[l + a,AyJ, (18) 
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Zy^ = 4°^[l + a,AyJ, (19) 



V4 



with 



Za. = Z'^l [1 + asAA^] , (20) 
Za, = Z^^lll + asAA,], (21) 



UQ = l-Zg^ I .^^nik, A = 0) = 1 - 1.2976(4)a, (22) 





/•TT , 




L (; 




Av,. 


2 


2 


A^i 


Av.2 


2 
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Av. = ^ + ^ + A,^, (23) 
Aa, = ^ + ^ + A,^,„ (24) 



where the integration in (p2D is performed numerically using the Monte Carlo 
integration routine BASES 0. A^ is the difference at the one-loop level between 
the lattice and continuum on-shell wave-function renormalization factors, and Ap 
(r = Jfj,,jfj.j5) are a similar difference for the vertex functions. We remark that 
A^ and Ap are functions of Ei'\ Ef^ and r. In the following two sections we 
present our calculation of A^ and Ap, and examine their pole mass dependences. 

3 Quark self-energy 
3.1 Lattice results 

The one-loop diagram for the quark self-energy is shown in Fig. 1, where ex- 
ternal quarks have zero spatial momentum. Using the lattice Feynman rules in 
Sec. 2 we can write the one-loop contribution as 

^E(^)(p4,m„,r) = ^ /" -^-^J^(A;,p4,m„,r) (25) 

Uq Uq J — it \^T^ ) 

with 

I^{k,p4, niu, r) = AttCf Vp{pi + /c/2, r)S{pi + k, niu, r)vp{p4 + k/2, r)D{k, A), (26) 

p 

where Cp = 4/3 is the quadratic Casimir invariant for the fundamental represen- 
tation of SU(3) group. The one-loop self-energy T,^^^ consists of the kinetic and 
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mass parts; 

S(^)(p4,m„,r) = i74sin(p4)E^^)(p4, m„,r) + J:^^\pi,mu,r). (27) 

With this expression the inverse of the quark propagator up to the one-loop level 
is given by 

5'"^(P4, m„, r) = uo\ Z74sin(p4) 



1 z^p\p4,mu,' 

Uq 

+2rsin2 (^4/2) - ^^^(p^, m„, r) \ . 



(2^ 



Since the quark mass is additively renormalized at the one-loop level due to the 
chiral symmetry breaking term in the Wilson quark action, on-shell condition for 
massless quark takes the form S~'^{p4 = 0, ml, r) = 0. Here the critical quark mass 
measures the magnitude of the additive renormalization. It is determined by 
the integral equation 



S«(p4 = 0,m:,r) 



d^k ^ ^ [m^ + 2rAi] AeW +rA4 
-(27r)4 4Ai{A4 + K + 2rAi]'} 

where, following the notation of Ref.0, 

Ai =^sin2(A;j2), 

a 

A4 = ^sin2(A;„), 

a 

A6(r) = (l + r')Ai -4. 
Using we can transform (pHj) into the following form 



(29) 
(30) 



(31) 
(32) 
(33) 



S ^{P4, ruu, r) =Uo{ i74sin(p4) 



1 ~ -j^S\p^^^u,r) 

) + ^S«(P4 = 0,m^,r)| 



Uq 



+ mu-ml 



+2rsin2 (^4/2) - ^^^(p^, rriu, r 



(34) 



Up to the one-loop level this expression is equivalent to 



S (p4, m„, r) = Uo< i74sin(p4) 



1 - ^^'p\P4,^u, 
Uri 



+2rsm^ (P4/2) - ^S^^(p4,m«,r) \ , 



(35) 



where 



rhu = mu- ml, (36) 
E(^)(P4,m„,r) = S«(p4,m„,r) - E^(^)(0, 0, r), (37) 

with ml = S^)(0, 0, r). We should note that, defining m„ as the bare quark mass, 



(|35D satisfies the on-shell condition for massless quark up to the one-loop level. We 
consider that m„ corresponds to the quark mass non-perturbatively determined 
from vanishing pion mass in Monte Carlo simulations. In the following analysis 
we use rhu as the bare quark mass. 

The pole mass E is given by the pole of the quark propagator; ^^^(^4 = 
iE,mu,r) = 0. At the tree-level the equation 

S'\p^ = iE^^\mu,r) = Mo{-74sinh(E(°)) +m„ - 2rsinh2 (^^°V2)} = (38) 

determines 

= In 



rhu + r + Jmi + 2rm„ + 1 



(39) 



1 + r 

One-loop correction is obtained by finding the pole of the quark propagator (^51). 
Defining the one-loop term by 

^ = ^(o) + ^^(i)(^(o)^^)^ (40) 

we find 



E«(E(°),r) 

sinh(EW)S«(p4 = iE^°\mu,r) - SW(P4 = iE^^\mu, 



(41) 



cosh(E(o)) +rsinh(E(o)) 
For later convenience we also present the expression of the one-loop correction to 
the pole mass including the tadpole contribution, which is easily obtained from 
(ID with the aid of (|23), 



„u m m m^ sinhC^^ + 2rsinh2 (e(o)/2) 

^ ' ^ V . ; cosh(E{o)) +rsinh(E(o)) ° ' ^ ^ 
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where u^^^ = -1.2976(4) from (H). In Table |T| numerical values of E^^\E^^\r) 
evaluated with r = 1 using BASES are given for representative values of E^^^ . 
The numerical accuracy is better than 2%. To show the magnitude of the tadpole 
contribution we also present in Table |l] the values of E"^^\E^'^\r) obtained from 
(P). Fig. 2 illustrate the dependence of E^'^\ 

The on-shell wave-function renormalization factor is defined as the residue of 
the quark propagator; 



(43) 



In terms of the expression (|35|) is written up to the one-loop level as 



[E^ ', r) = Uq [cosh(£/) + rsmh{E)\ ^ 



where E is given by (^) and 



(44) 



P474=iE(0) 



S(^)(p4,m„,r) = i74sin(p4)S^^)(p4,m„,r) + S^)(p4,m^„r). (45) 

Even at the tree-level the on-shell wave-function renormalization factor suffers 
from finite m^a corrections as is clear from 

4°^~'(E(o),r) = Mo [cosh(E(o)) +rsinh(E(°))] . (46) 

Factorizing the tree level contribution we obtain the following expression 

ZAE^'\r) ^ a. sinh(EW)+rcosh(EW) 
4°^"'(E(o),r) M^cosh(E(o))+rsinh(E(o)) ^ ' ^ ^ 

as 1 9S(^)(p4,m„,r) 



P474=« 



cosh(E(o)) + rsinh(E(o)) i9(p474) 
Applying the power counting rule we note that the right-hand side has an IR 
divergence for A ^ arising from the third term. To extract the terms independent 
of A we consider subtracting from the integrand of the third term an analytically 
integrable expression which has the same IR behavior in the region of small loop 
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momentum k; 



cosh(E(o)) +rsinh(E(o)) ^^(^474) 
d'^k dl4k,p4,E^^^) 

d^k d 



P474=ii?('') 



+ 



(27r)4 Z9(P474 

where the IR behavior of 



P474=ii?('') 

/^(fc,p4,m„,r) f fi. „ f{o)n 

cosh(E(o))+rsinh(E(o)) ) 



is the same as that of 
d 



p474=i_B(0),A=0 



(48) 



P474=iE(0),A=0 



(49) 



P474=iE(0),A=C 



(50) 



d I^{k,p4,mu,r) 

id{pi^i) cosh(E(o)) + rsinh(E(o)) 

In this expression the IR divergence is transferee! to the first term, and in conse- 
quence the second term is finite. For the candidate of we try the continuum 
counterpart of the integrand J^/[cosh(i?'^'')) +rsinh(i?*^°^)] replacing Vp{pi + k/2^ r) 
by Vp, S{pi + k, r) by S{pi + k, E^o)) and D{k, A) by D{k, A) in (0), 

i4k,p,,E^'>'^) = 9{A' - e)4nCFY.^pS{P^ + k, E^^^)vpD{k, X), (51) 

p 

where the domain of integration is restricted to a hyper-sphere of radius A, not 
exceeding vr, for convenience of an analytical integration. It is apparent that in 
the limit of a -> (||) and (|50|) have the same IR behavior. At a finite lattice 
spacing, however, the IR behaviors of the two integrands are different due to finite 
rUga corrections. 

Let us examine the IR behaviors of the denominators of (^) and (0), where 
we transfer the Dirac structure in the denominators of quark propagators to the 
numerators. For the continuum case we find 



D-\k, A = 0) [s-\p^ + k, E(°))5-i(-(p4 + k), E(°)) 
2 



k' 



(52) 
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with the use of the on-shell condition = iE'^^\ For the lattice case we expand 
the quark and gluon propagators around A; = 0, obtaining 

cosh(E(°)) + rsinh(E(°))] D-^{k, A = 0) 
X {S~^{.P4. + A;,m„,r)S'"^(-(p4 + A;),m„,r)} 



cosh(E(°)) + rsinh(E(°))] + O^k" 
X {i2A;4sinh(E(°)) [cosh(E(°)) + rsinh(E(°))] + k'^(l + rsinh(^(°))) 
+A;2sinh(E(°)) [(2 - r2)sinh(E(°)) + r(cosh(^(°)) - 1)] + 0(P)}^ , (53) 

where the on-shell condition of external quark is used. While the two expressions 
above show a different IR behavior, the difference can be absorbed, aside from an 
overall factor, by replacing the pole mass in (0) by 

cosh(E(°)) + rsinh(E(°)) 



m 



sinh(E(°))- 



(54) 



l + rsinh(E(o)) 

Here, we do not take account of the last term in ( |53D because it does not give 
leading order contributions either for the case of A;4 = or of A;4 7^ for small 
k. It is straightforward to check that the remaining overall irtqa corrections of the 
denominator are precisely canceled with those arising from the numerator with 
A; = 0. 

Consequently we take for the integrand the following expression, 

i^{k,pi,m) = e{ls? - k'^)A'nCF^VpS{pA + k,m)vpD{k, A). 



(55) 



A simple calculation gives 

f"^ d^k dl^{k,p4,m) 



. (27r)4 ^9(P474) 



47r V 
-6 In 



-2 In 



A^ 



A2 



3A^ 



P474=«m 

9A 



4m^ 2711? 



3A 



IS? + Am? H -(A^ + Am 



A + ./A2 + 4m2 



2m 



whose massless limit is 

r'^ d^k dl^{k,p4,m) 



lim 

m— >0 



(27r)4 z9(p474) 



P474=4m 



Cf 
An 



-In 



A^ 



- 2 In 



A^ 



(56) 



- o ^57) 
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where A is assumed to be less than m. 
3.2 Continuum results 

We turn to the calculation of the on-shell wave-function renormalization factor 
using the continuum NDR scheme instead of the lattice regularization scheme. For 
the one-loop contribution to the quark self-energy shown in Fig. 1 the continuum 
Feynman rules in Sec. 2 give the expression 



POO ~ ~ 

asTy\p,m) ^ as / 47rCF A:, m)i;p.D(A:, A), 

J-oo [Z-K) p 



(58) 



where D is the reduced space-time dimension which is parameterized by e as 

D = 4 - e, e > 0. (59) 

Since this dimensional reduction procedure prevents us from taking zero spatial 
momentum for the external quark state before the loop integration, we perform 
the calculation of the on-shell wave-function renormalization factor in a Euclidean 
invariant way. 

Up to one-loop corrections the quark propagator takes the form 

S~ {p, m) = if) + m — as S(^)(p,m). (60) 

In terms of this expression the on-shell wave-function renormalization factor Z^^ 
is written up to the one-loop level as 

aEW(p, m) 



Z^{m) = l-as 



idj) 



(61) 



At the tree level we find 



(62) 



Performing the integration in an elementary way we obtain the one-loop correction 



idj) 



9l 

An 



(^--'j + ln - In 



In 




-2 In 


A2 


-4 













(63) 
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where the pole term (2/e — 7 + ln|47r|) should be eliminated in the MS scheme. 
We note that this result is independent of the choice of spatial momenta for the 
external quark state because of Euclidean invariance in the continuum theory. 

3.3 Relation between continuum and lattice wave-function 
renormalization factors 

Using the results for the lattice and continuum wave-function renormalization 
factors obtained in the previous subsections, let us find the correction factor which 
connects the two factors. From (^) and (^) the tree-level relation is 

= % [cosh(EW) + rsinh(E(0))] zf\E^'\r), (64) 
where we take m = E^^\ Up to the one-loop level we obtain the following expres- 



sion 



Z^°°t(m) = Mo [cosh(E(°)) + rsinh(E(°))] [l + a,A^(E(°), r)] Zlj^'\E^^\ r) , (65) 
where, from ( ^Tj) and (|6lD, 

sinh(EW)+rcosh(igW) (,)iatt 
cosh(E(o)) +rsinh(E(o)) ^ ' ' 



-1171 



1 9S(^)^''"(p4,m„,r) 



, (66) 



cosh(E(o)) +rsinh(E(o)) ^^(^474) 

with m = E^^\ From (^2]), ([65| ) and (|66D the tadpole contribution to the lattice 
wave-function renormalization factors is explicitly expressed as 

A;(i?(°\r) 

= A^(EW,r)+4^) 

[sinh(E(°)) + rcosh(E(°))] [sinh(E(o)) + 2rsinh2 /2)] 



[cosh(E(o)) + rsinh(E(o))]^ 



(67) 



where u^^^ = —1.2976(4) from (P^. For the coupling constant we can take either 
Q,cont Os^'^^^/uq because the difference is of order a^. From (0) and (^) we 
observe that the IR singular terms in (pBI) for A ^ are precisely canceled. We 
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also note that the mass singularities occuring at E^^^ — > in individual terms of 
(|66D also cancel, which assures us that A^(E(o) , r) is finite even in the massless 
limit. 

In Table |l] we present numerical values of A^(E^^\r) evaluated using BASES 
with an inaccuracy of less than 2% for representative values of the pole mass 
for the r = 1 case. The values of A'^{E^^\r) defined in ( |67D are also given in 
Table to demonstrate the magnitude of the tadpole contribution. Our result for 

evaluated at E^^^ = with r = 1 is 

A;(E(°) = 0, r = 1) = A^(E(°) = 0, r = 1) + = -1-37, (68) 



which is consistent with Asi/(37r) for r = 1 in Ref. ||T0| within the error in the 
numerical integration. The pole mass dependence of A^ is shown in Fig. 3. We 
observe that the m^a corrections give large contributions in the region 

4 Vertex functions 
4.1 Lattice results 

At the tree level the lattice vertex functions for F = 7^, 7^75 (/x = 1, ■ ■ ■ , 4) do 
not suffer from rriqa corrections, which contrasts with the wave-function case. Up 
to the one-loop level the vertex functions are written in the following form 

Ar{Er\Er\r) = r + ^,A?{Ei'\Ei'\r) (69) 

where e[^^ and E^f^ are the pole masses for the external quarks. The one-loop 
vertex corrections Ap'' for F = 74,7475 and F = 71,7475 {i = 1,2,3) are obtained 
by evaluating the diagrams shown in Fig. 4(a) and Fig. 4(b), respectively, on 
condition that the external quark and anti-quark are on-shell with zero spatial 
momentum. 

We first consider the case of F = 74, 7^75 for which external states are quarks. 
Using the lattice Feynman rules in Sec. 2, the amplitude corresponding to Fig. 4(a) 
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is expressed as 



a. 



J-n (27r)' 



lAk,E?\E^^\r) 



with 



(70) 



Ir{k,E?\Ei'\r) = 4nCFj2^piP'4 + k/2,r)S{p'^ + k,m2u,r) 

p 

xTS{p4 + k, mi„, r)vp{pi + k/2, r)D{k, A), (71) 

where and i?2°^ are expressed with rhiu and m2u, respectively, as in (p9|). The 
vertex correction ([7^) has IR divergences for A ^ as can be seen by the power 
counting. Expanding (|7l|) around k = we extract the IR singular part: 

I,,ik,E['\Ei'\r) 

I [AnCFD^k, \)S{p4 + k,rhiu,r)S{-{p4 + k),miu,r) 



xS{p^ + k,m2u,r)S{-{p^ + k),m2u,r)] 

1 + 



74 



2sinh(E{°^)sinh(E, 



(Oh 



cosh(E;°)+4°^) 



cosh(Ef ^ - 4°^) - rsinh(Ei"^ + E^' 



(0) , 7^(0) ^ 



+ 



2sinh(Ef^)sinh(E^' 
l-r2 



(0)^ 



1 + 



(0) , 7^.(0) ^ 



cosh(El"^ + E'2 



cosh(Er - ^^'0 - rsinh(E{°) + 4°^^ 



2 

+0(fc), 

I,^,,ik,E?\E!t\r) 

I ['i7rCFD{k,\)S{pi + k,7hiu,r)S{-{p4 + k),miu,r) 
xS'(p4 + k,rh2u,r)S{-{p'^ + A;),m2„,r)] 
7i75 2sinh(E{°^)sinh(£;f^ 



1 _ l±z:!cosh(E;°) + 



(0)^ 



l-r2 



cosh(E{°) - 4°^) - rsinh(Ei"^ + E)- 



(0) , 77.(0). 



+7^7574 



2sinh(Er)sinh(Er) < 1 



-cosh(£;r + ^r) 



(72) 
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1-r' 



+Oik), 



(73) 



where, from 



D ^{k,X)S ^{pi + k,miu,r)S ^{-{p^ + k),miu,r) 

fc2 + + 0{k^) 

{22fc4sinh(Ef ^) [cosh(E{°^) + rsinh(Ef^; 
+A;2(l + rsinh(E{°))) + 0(A;2)} 
X {i2fc4sinh(Ef ^) [cosh(Ef ^) + rsinh(E^°^; 
+A;2(l + rsinh(Ef )) + 0(A;2)}, 



X 



(74) 



with the use of the on-shell conditions = iEf'' and p'^ = iElfK In order 
to extract the terms independent of A we design to subtract from the integrand 
Ir{k, Ei'\ E2'\r) an analytically integrable expression Jp which has the same IR 
behavior near k = 0; 

d^k ~ 



n (27r)4" 
- d^k 



(75) 



(27r) 



h{k,Ef\Ef\r) - h{k,fn^,m2] 



A=0 



where the first term in the right hand side contains the IR divergence and the 
second one is IR finite. For the candidate of Jr we take 



ir{k,mi,m2) = 9{A^ - k'^)A7rCFY.^P^iP'4 + ^^^2] 

p 

xrS'(p4 + k, fhijVpDik, A), 



(76) 



with 



mi,2 = sinh(Er2' 



(o)^cosh(Eg)+rsinh(Ei°j) 



1 + rsinh(Ef2^) 



(77) 
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as we did for the case of the on-shell wave-function renormahzation factor in Sec. 3. 
The IR behavior of Jr is shown by expanding (^) around k = 0, 

I ^^{k, nil, 7X12)1 4:TcCFD{k, X)S{p4 + k,mi)S{-{pi + k),mi) 

xS{j)'^ + k,m2)S{-{j)'^ + k),m2) 
= 74 [4mim2] + [-Smima] + 0{k), (78) 

/^^^5(/c,mi,m2)/ A-nCpD^k, X)S{pa + k,mi)S{-{pA + k),mi) 

xS{p'^ + k,m2)S{-{p'^ + k),m2) 
= 7i75[-4mim2] + 0(/c), (79) 



where, from (E2\ 



D-\k, \)S-\p^ + k, mi)S-\-{pi + k),mi] 
xS'\p^ + k, Th2)S~\-{p^ + k),Th2) 
k'^ + A^l \i2k4,mi + \i2k4m2 + k'^^ 



^0) 



with the use of the on-shell conditions p4 = irhi and p'^ = 17712- 

Comparing (|72D with (^) we find that each term proportional to 74 and 1 in /^^ 
shows different IR behavior from that in J^^ at the finite lattice spacing a ^ 0. This 
is not due to the operator mixing between the local vector current and the local 
scalar density but ascribed to the lattice cut-off effects, because it is known that 
the local vector and axial vector currents do not mix with other dimension three 



operators for the Wilson quark action|jTO|. (In this paper we do not consider the 
contributions of the higher dimensional operators to the lattice vertex functions 
(|69D as the operator mixing.) Here it is noted that following our definition of the 
renormahzation constant in (|1^) we should combine the term proportional to 74 
with that proportional to 1 in /^^ and I^^ using the equations of motion for the 
external quarks on the lattice 



(i74sin(j94) + TTiiu + 2rsin^ (P4/2)) u{p4) = 0, 
u{p'a) (^74sin(p'4) + 7h2u + 2rsin^ (p'4/2)) = 0, 



(81) 
(82) 
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and in the continuum 



{iim + mi) u{pi) = 0, (83) 
(^74p1 + rh2) = 0, (84) 

where ^(^4) is the Dirac spinor for the incoming quark state and ulp^) that for 
the outgoing quark state. With this procedure we find that the IR behaviors of 
and are same, which is expected from a point of view that the IR behavior 
should be independent of the ultra-violet regularization scheme. This is the case 
that we can match the vertex correction for F = 74 on the lattice to that in the 
continuum irrespective of the IR regularization scheme. For F = 7^75 the combined 
contribution of the terms proportional to 7^75 and 7^7574 in J^.^g of ( |73D with the 
aid of the equations of motion for the external quarks (^) and ( p2D has the same 
IR behavior with J^.^g 

of (©. 

In the expression (^) we have carried out the first integral analytically and the 
second one numerically using BASES. The analytical integrations of Ir{k, mi, 7712) 
for F = 74, 7j75 give the following results 
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m2 



mi — 7712 
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A2 + 4mf A 
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A 
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4mfm2 
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m2 
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and 



(27r)4 ^^^^ 
A2 



(/c,mi,m2)/ 



47r 



7i75 
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m2 



A2 



+ 



^TTii 71717712 J 2 y2mfm2 
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71717712 J 12 



mi — m2 
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' A2 + 4m? A / 

-4^^ : + ^JK" + 4mf 
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— (A^ + 4m?^ 
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6 In 



2mf 



A + a/A2 + 4mf 



2mi 



+ (mi 



m2) 



(86) 



where we use the equations of motion (^3|) and (Q) for the second expression in 
(|S5D. We remark that above expressions do not diverge in the hmit of m2 mi. 
The resuhs for mi = m2 are presented in Appendix. 

Let us turn to the second case F = 7j, 7475 for which a quark and an anti-quark 
have to be chosen for the external states for a non- vanishing matrix element. From 
Fig. 4(b) we obtain (|70|) with, however, 

h{KEf\Ef\r) = ^'KCFY.^p{-p'^ + k/2,r)S{-p'^ + k,77i2u,r) 

p 

XF^(P4 + k, mi„, r)vp{pi + k/2, r)D{k, A). (87) 



Around k = the integrand (|87D behaves as 



I [4:nCFD{k,\)S{p4 + k, 777iu,r)S{-{pi + fc),mi„,r) 



xS'(-p4 + k,rh2u,r)S{-{-p'^ + k),rh2u,r)] 

1 + 



+ 



2sinh(E{°^)sinh(£;^°^) |l 



cosh(Ef) 



(0)^ 



cosh(E}^^ - E^^O + rsinh(E{°^ + ^2 
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where 



+7^74 



(0)^ 



-1 + 



1 



cosh(£;i°^ - E^"') + rsmh{Er + E^) 



1 + 



(0) 



cosh(£;f ^ + E^; 



2 

+0(fc), 



.(0)^ 



A475 -^1°'' 1 -^2°'' ) ^) 

/ [47rCF-D(fc, A)S'(p4 + fc, rhiu, r)S{-{p4 + k),miu, r) 
X'S'(-P4 + k,m2u,r)S{-{-p'4^ + k),m2u,r)] 

2sinh{El%mhiEi'^) (-3 + l±i:!cosh(£;f ^ + E^'^) 



7475 



1 -r^ 



cosh(£;i°^ - E^) + rsmh{Er + E'^ 



(0) 



.(0)^ 



+75 



2sinh(£;i°^)sinh(£;r) { -3 



(0)^ 



+0{k), 



cosh(£;i°) - - rsinhiEl"" + E^^ 



.(0) 



1 + 



(0) 



coshiEf"^ + E^">) 



(Oh 



(89) 



D ^{k,X)S ^{p4 + k,miu,r)S ^{-{p4 + k),miu,r) 
x5'"^(-P4 + k,m2u,r)S~'^{-{-p'4 + k),m2u,r) 
= [A;' + A' + 0(A;^)" 

X {i2A;4sinh(Ef^) [cosh(Ef ^) + rsinh(Ef 

+A;^(l + rsinh(£;f )) + 0(A;4^)} 

X {-i2A;4sinh(4°^) [cosh(4°^) + rsinh(Ef 

+A;^(1 + rsinh(£;f )) + 0(A;|)} , (90) 

with the use of the on-shell conditions p4 — iE[^^ and P4 = i£^2°'* ■ 

We introduce a counter term Ir{k,rhi,rh2) to regularize the IR singularity of 
the integrand /p in the vertex correction, 

ir{k, 1711,1712) = ^(A^ - fc2)47rC^^{)p^(-p; + /c,m2) 
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xTS{p4 + k, mi)vpD{k, A) 



(91) 



with mi 2 defined in (|77D. The IR behavior of Jr is 



I ^.{k,fhi, 7X12)1 4:-7TCFD{k, X)S{p4 + k,mi)S{-{pA + k),rfii) 
xS{-p'^ + k,fh2)S{-{-p'^ + k),m2) 
7i [4mim2] +0{k), 



(92) 



'■7475 V 



(fc, mi, m2)/ AnCpDik, \)S{pi + k, mi)S'(— (^4 + A;), mi) 



xS'(-p4 + k,m2)S{-{-p'^ + k),m2) 
7475 [-4mim2] + 75 [-8mim2] + 0(A;) 



(93) 



where 



^-^(A;, A)5-^(p4 + k, mi)S-\-{pi + k), 
xS-\-p'^ + k,m2)S'\-i-p'i + k), ma) 



i2kimi + /c^ —i2k4m2 + /c 



(94) 



with the use of the on-shell conditions ^4 = imi and p'^ = im2. 

As we did in the case of the vertex corrections for F = 74,7175, we make a 
one- loop matching of the vertex corrections for F = 7^, 7475 on the lattice to those 
in the continuum with the aid of the equations of motion for the external quark 
and anti-quark states on the lattice 



(i74sin(p4) + mi„ + 2rsin^ (^4/2)) ^(^4) = 0, 
v{p'^) (-i74sin(p'4) + m2„ + 2rsin^ (p'4/2)) = 0, 



(95) 
(96) 



and in the continuum 



{i'^iPi + mi)u{pi) -- 

vip'4) (-^74^4 + "i2) 



0, 



(97) 
(98) 



where u{p4) is the Dirac spinor for the incoming quark state and v{p'^) is that for 
the incoming anti-quark state. 



23 



Using Ir{k, mi,m2) of ( pTD we can decompose the vertex correction A^\Ei'\ E!f\r) 
into an IR divergent part and a finite one as in ([75|) . The IR divergence residing 
in are calculated analytically. The finite term is evaluated with numerical 
integration with the aid of BASES. The results of an analytical integration of (|9lD 
for r = 7j, 7475 are as follows, 
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where we use the equations of motion (^3) (|98D for the second expression in 
(|100|) . The expressions in the hmit of 7/7.2 —>■ ^1 are hsted in Appendix. 

4.2 Continuum results 

We repeat the calculations of the vertex functions for T = 7^,7^75 (/x = 
1, • • • , 4) in the continuum NDR scheme. The vertex functions up to the one-loop 
level are written as 



Ar(mi,m2) = T + asA^jl\mi,m2). 



(101) 



The one-loop contributions Ap-* for T = 74,7^75 and T = 7^,7475 {i = 1,2,3) 
are represented in Fig. 4(a) and Fig. 4(b), respectively. We perform the loop 
integrations imposing the on-shell conditions on the external quark and anti-quark 
states in a Euclidean invariant way, and after the integrations set the spatial 
momenta of the external quark and anti-quark equal to zero. 

For the case of F = 74,7^75 the continuum Feynman rules in Sec. 2 give the 
following expressions for Fig. 4(a), 



asAp^(mi,m2) = a 



00 (j^Dj^ ^ 

s I -j—-^A'kCf VpS {p' + k, 1712) 

00 



(27r)^ 

xTS{p+ k,nii)vi,D{k,X), 



(102) 



where D is the space-time dimension which is reduced from four by e to regularize 



ultraviolet divergences and F is defined as 7^ or 7^75 for /i = 1,2, ■ ■ ■ , D. The 



25 



on-shell conditions for the external quarks are written as 

[i-p + mi)u{p) = 0, 
u{p){i-^'+ 777,2) = 0, 



(103) 
(104) 



where u{p) is the Dirac spinor for incoming quark state and 'u(p') is one for outgoing 
quark state. After carrying out the integration in (|102|) we take 7^ = 74, 7^75 = 
7i75! P = (0, 0, 0, irrii) and p' = (0, 0, 0, 17712). The final results are given by 
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(106) 



where the pole term (2/e — 7 + ln|47r|) should be eliminated in the MS scheme. 
The second expression in ( |105| ) is obtained using the equations of motion (|103|) 
and ( |104|) . The expressions of (|105|) and (|106|) for the case of 7772 mi are given 
in Appendix. 

Another case is F = 74,7475. The one-loop diagram shown in Fig. 4(b) is 
written as 



asA[!'' (7771, 7772) 
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xTS{p + k, mi)vpD{k, A), 
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(107) 



where the external quark and anti-quark are on-shell 



{i-jf) + mi)u{p) = 0, 
v{p'){-if + 1712) = 0, 



(108) 
(109) 



with v{p') the Dirac spinor for outgoing anti-quark state. Performing the integra- 
tion we obtain 
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where the pole term (2/e — 7 -|- ln|47r|) should be eliminated in the MS scheme. 
We use the equations of motion ( |108| ) and ( |109| ) for the second expression in ( [1111 ). 



4.3 Relation between continuum and lattice vertex func- 
tions 



From (^) and ( |101|) we obtain the relation between the lattice vertex correc- 
tions and the continuum ones up to the one-loop level 



(0) r^iO) 
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with 

Ar{E['\E!f\ r) = Ai^')™°Vi, ^2) - Af^'^'^^V^, , r), (113) 

where we take mi = e[^^ and m2 = E!f\ In Appendix we hst the expression for 
A[!^™''* and A't^^^^^ in the hmit e[°^ = E!f^ ^ or ^ 0. They show that 
Ap and Ap have the same singular structures for X ^ 0, E^ = E2 ^ 
or Ei'^ — > where we assume A < Ei^\E2'\ Thus Ar in ( |113| ) is IR finite, and 
also finite in the limit of e[°^ = E^^^ ^ or e[^^ 0. 

Numerical values of Ar{Ef'\ E2^\r) for F = 7^ evaluated using BASES with 
an accuracy of better than 2% are tabulated in Table |^ for representative values 
of the pole masses Ei'^ and E^f^ for the r = 1 case. Our results for A^. and A^^ 

(0) _ p(0) 



at El = E2 — are consistent with A^^/ (Svr) for r = 1 in Ref. [jTO[ within the 
error in the numerical integration. Fig. 5(a) shows the i?2°'* dependence of A^^ and 
A^4 for the case of Ei'^ = -^2°'', and Fig. 5(b) is the same as Fig. 5(a) for Ei'^ = 0. 
For both cases we observe that A^. and A^^ are roughly consistent in the region 
-£^2°^ ^0-01) while as E^^ becomes larger the absolute value for A^^ increase and 
that for A^^ decrease, indicating large m^a corrections. 

Numerical data of At{e[^\ E^\r) for F = 7^75 are listed in Table |^. The 
values of A^.^^ and A^^^^ at Ei'^ = E^f^ = show an agreement with A^^^5/(37r) 
for r = 1 in Ref. within the error in the numerical integration. Our results for 
A^.^5 and A^^^^ are plotted in Fig. 6(a) for the case of e[°^ = E^'^ and in Fig. 6(b) 
for Ei'^ = 0. We observe that A^.^^ and A^^^^ show an E^'^ dependence similar 
to that for A^^ and A^., respectively. 

5 Renormalization factors for vector and axial 
vector currents 

Since we have completed the calculations of A^ and Ar (F = 7/^,7^75) we 
are now ready to discuss the magnitude of finite rriqa corrections in the one-loop 
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contributions of renormalization factors for vector and axial vector currents. Com- 
bining the results for and Ar we obtain values of Ay^ and A^^ defined in ( ^31) 
and (H). 

We show the results for Av^X-^i^"*) -^2°\ Fig- 7(a) for Ef'' = E^^ and in 
Fig. 7(b) for Ef''^ = for the r = 1 case. For both cases the values for Ay^ are 
close to those for Ay^ below E2^<^0-01, beyond which Ay^ decreases in magnitude, 
whereas Ay. increases. In the heavy quark mass region ^0(1-2), where 
current 6-quark simulations are performed, the values of Ay and Ay^ deviate by 
about 100% from those at ^ = 0. 

For the axial vector current the results for Aai,{e[^\ E!f\r) are plotted in 
Fig. 8(a) for e[^^ = E^f^ and in Fig. 8(b) for Ei'^ = 0, again for the r = 1 case. 
We observe that the characteristic features for the E!f^ dependence of A^^ are 
similar to those of Ay^, where A^^ corresponds to Ay^ and A^^ to Ay. We note 
that A^. and A^^ also suffer from a 100% m^a correction in the heavy quark region 
^0(1-2). 

6 Large quark mass limit of renormalization fac- 
tor for heavy-light axial vector current 

It is instructive to examine the quark mass dependence of the renormaliza- 
tion factors in the heavy quark region toward the static limit comparing with 
the previous results for the static effective theory [[11] , |l^ and the nonrelativistic 
QCD(NRQCD) UlSj [T^ . In the static case there exists only the calculation of the 



renormalization factor for the fourth-component of the static-light axial vector 



current using a static heavy quark and a massless Wilson quark ITT], |T2[. Thus in 
this section we consider the renormalization factor of the heavy-light axial vector 
current. It is expected that our Wilson results should agree with the static ones 
in the large quark mass limit. 

In Table ^ we present numerical values for E^^^^ {E2^\r), A'^(E2^\r) — 
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2/7r In 



E. 



(0) 



and A,y^^^{Ei'^ = 0,E2'',r) evaluated with r = 1 for representative 



n(0) 



.(0) 



values of E2^\ where the divergent part for the heavy quark mass limit E2' — > 00 
in is subtracted as before. For the values at E^f^ = 00 we quote the results in 



Ref. Since the static results for the one- loop correction to the pole mass and 



that to the wave-function renormalization factor in Ref. |11] contain the tadpole 



contributions, we present our results for £'2'^^'' defined in p^ ) and AL in (|67|). 



latt 



Figures 9, 10 and 11 show the E^^^ dependence of E2"'''"~^\ In 
and A^4^g respectively. For the comparison we also plot the results for the NRQCD ||T3| 
where we use the following correspondences between our notations and theirs: 



/(I) 



latt 



latt 



a; - 2/n In 



E. 



(0) 



^7475 ' — ' 



-AixaA 



47r4^", 



(114) 
(115) 
(116) 



with £'2°'* the continuum bare quark mass in common between our Wilson results 
and the NRQCD ones. We observe that our results become roughly consistent with 
the static limits around E'2^ 5, while the NRQCD results become closer to the 
static limits rather slowly having larger values in magnitude than our results. It is 
a reasonable feature that the NRQCD results deviate from our Wilson ones toward 
the lighter quark masses, since the NRQCD is well-defined only in the heavy quark 



mass region. In terms of AJ^(_E'j°'' = 0,r) = —1.37 in (|68l), A^,,(-5'^°'', r)— 2/?? In 



and A^4^5(eJ°^ = 0,Er,r) we construct A^^(Er = 0,E^''^r) - l/vr In 
following (^), whose E^^"* dependence is given in Fig. 12. Both results for the 
Wilson quark action and the NRQCD show a smooth approach to the static limit, 
while the latter has larger absolute values than the former over Ef^^^l. 



.(0) 



.(0) 



.(0) 



4°) 

Ef^ 



7 Conclusions and discussion 



In this paper we have calculated the one-loop contributions for the renormal- 
ization factors of the vector and axial vector currents including finite quark mass 
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corrections using the Wilson quark action. We have demonstrated that the one- 
loop contributions suffer from a very large correction of 0(100%) for the heavy 
quark masses of order unity in lattice units, which corresponds to the situation 
of current 6-quark simulations. This fact tells us that the m^a corrections should 
be incorporated even at the one-loop level for the renormalization factors of the 
weak operators containing the 6-quarks in current numerical simulations using the 
Wilson quark action. We have also checked that the one-loop contributions for the 
renormalization factor of the heavy-light axial vector current show an agreement 
with the static results toward the heavy quark mass limit. 

In this work our investigation has been concentrated on the Wilson quark 
action, which is the most naive, unimproved, quark action on the lattice. We 
wonder how large the m^a corrections are for improved quark actions at the one- 
loop level. This point should be examined through a calculation similar to ours 
including finite quark mass contributions, which we leave for future investigations. 
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Appendix 

In this appendix we present expressions of the continuum one-loop vertex cor- 
rections and the integrals of the counter terms introduced to regularize the IR 
divergence of the one-loop lattice vertex corrections for some specified cases of the 
pole masses. 

We first list expressions of the continuum vertex corrections Ap"* (mi, 777.2) (r = 
7^, 7y^75) for the case of m2 — > 777.1 and r77i — > with the assumption A < r77i, 7772- 
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From ( |105| ), (|106|) , ( |110| ) and ( |1 1 1| ) we find that the vertex corrections are written 
in the following forms 
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where the second expressions in ( 118 ) and (|122|) are obtained using the equations 
of motion for the external quarks ( p.03| ) and (|104|) , and those in ( p.2CI| ) and (|124|) 
are with the equations of motion for the external quark and anti-quark (|108|) and 

We also enumerate expressions of the integrals of the counter terms to the 
lattice vertex corrections J^^ d^k/ {2n)^Ir{k, 1711,1712) (T = 7^, 7;i75) for the case of 
m2 fhi and mi with the assumption A < mi,m2. From 
and ( |10(J| ) we obtain the following expressions 
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where the second expressions in ( |126| ), ( |13CI| ) and ( |134| ) are obtained using the 
equations of motion for the external quarks (p3D and 



, and those in (|128D, 



( |132| ) and ( |136| ) are with the equations of motion for the external quark and anti- 
quark (pTf ) and 
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Tables 



Table 1: dependence of E^^') \E^^\ r) , (EW,r), A^(E(o),r) and 

A'^(E^^\r). The Wilson parameter r is chosen to be one. For and 

the error in the numerical integration is less than 2%. The values for and 
A^ are evaluated using (^) and ([671). 
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Table 2: (a) A^.{Ei ' , E2 ,r) and (b) A^^{E[ ' , E2 ,r) for combinations of El ' and E2 . The Wilson parameter r is chosen to 
be one. The error in the numerical integration is less than 2%. 
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Table 3: (a) A^.^^{E[ ' , E2 ,r) and (b) Ay^^^{E[ ' , E2 ,r) for combinations of E[ ' and E^ The Wilson parameter r is chosen 
to be one. The error in the numerical integration is less than 2%. 
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-1 


-3.25x10" 


-1 


-3.10x10" 


-1 


-3.13x10" 


-1 






1 


-4.83x10- 


-1 


-4.74x10- 


-1 


-4.89x10" 


-1 


-5.15x10" 


-1 


-5.21x10" 


-1 


-4.99x10" 


-1 


-5.46x10-1 




2 


-8.16x10- 


-1 


-8.20x10- 


-1 


-8.38x10" 


-1 


-9.27x10" 


-1 


-1.00 




-1.12 




-1.20 


-1.61 



Table 4: dependence of E2'^'>"^^\E. 



(1) 



latt , 



n{0) 



(0) 



- 2/71 In 



A7475 (-^1 



(0) 



0,E, 



(0) 



E. 



(0) 



and 



The Wilson parameter r is chosen to be one. The er- 



ror in the numerical integration is less than 



70. 



Ef^ 



1.0 
1.4 
1.8 
2.2 
2.6 
3.0 
4.0 
5.0 
00 



11 



TatF 



1.33 
1.53 
1.70 
1.81 
1.90 
1.96 
2.07 
2.13 
2.117 



a; - 2/7r In 



-7.75x10^^ 
-8.26x10^^ 
-8.52x10 
-8.81x10 
-8.76x10 
-8.82x10 
-8.85x10 
-9.12x10 
-9.05x10 



eP 



A, 



7475 



-4.83x10" 

-6.21x10" 

-7.55x10" 

-8.74x10- 

-9.78x10" 

-1.06 

-1.17 

-1.20 

-1.239 
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Figure Captions 

Fig. 1 One-loop diagram for the quark self-energy, k is the loop momentum and p 
is the external quark momentum. 

Fig. 2 One- loop correction to the pole mass on the lattice as a function of E^*^^ . The 
Wilson parameter r is chosen to be one. Open symbol denotes the value at 

= 0. 

Fig. 3 One-loop coefficient of the relation between the on-shell wave-function renor- 
malization factors on the lattice and in the continuum with NDR scheme as a 
function of E^^^ . The Wilson parameter r is chosen to be one. Open symbol 
denotes the value at E^°^ — 0. 

Fig. 4 One-loop diagrams for the vertex corrections, k is the loop momentum and p 
is the incoming quark momentum, p' denotes the outgoing quark momentum 
for (a) and the incoming anti-quark momentum for (b). 

Fig. 5 One-loop coefficients of the relation between the vertex corrections for 7^ on 
the lattice and in the continuum with NDR scheme, (a) is for the case of 
e[°^ = E^°^ and (b) for e[^^ — 0. The Wilson parameter r is chosen to be 
one. Open symbols denote the values at E^'^^ = 0. 

Fig. 6 One-loop coefficients of the relation between the vertex corrections for 7^75 
on the lattice and in the continuum with NDR scheme, (a) is for the case of 
e[°^ — E^^^ and (b) for e[°^ — 0. The Wilson parameter r is chosen to be 
one. Open symbols denote the values at E^'^^ = 0. 

Fig. 7 One-loop coefficients of the renormalization factors for vector currents, (a) 
is for the case of Ei^^ — £^2°'* and (b) for Ei^^ — 0. The Wilson parameter r 
is chosen to be one. Open symbols denote the values at E^^^ — 0. 
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Fig. 8 One-loop coefficients of the renormalization factors for axial vector currents. 



(a) is for the case of Ef^ = E2'' and (b) for E^'' = 0. The Wilson parameter 
r is chosen to be one. Open symbols denote the values at E^^^ = 0. 

Fig. 9 E2 dependence oi E2 (£'2 , t) for the Wilson quark action and the 
NRQCD[p^. The Wilson parameter r is chosen to be one. Dotted line 



.(0) 



(0) 



denotes the static result in Ref. 11 



Fig. 10 dependence of A'JE^^^r) - 2/7r In 



(0) 



E. 



(0) 



for the Wilson quark action 



and the NRQCD||13||. The Wilson parameter r is chosen to be one. Dotted 



line denotes the static result in Ref. 11 



Fig. 11 Elf^ dependence of A^^^^{Ei'^ = 0,E2'\r) for the Wilson quark action and 
the NRQCD||13||. The Wilson parameter r is chosen to be one. Dotted line 



denotes the static result in Ref. 11 



Fig. 12 ^ dependence of A'^^iEl"" = 0, ^) - V^r In 

action and the NRQCD||13||. The Wilson parameter r is chosen to be one. 
Dotted line denotes the static result in Ref. |]ll] . 



.(0) 
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E. 
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